Minimum energy needed to perform a quantum logical gate 



Julio Gea-Banacloche|2] 

Department of Physics, University of Arkansas, Fayetteville, AR 12101 
(Dated: February 1, 2008) 

A lower bound on the amount of energy needed to carry out an elementary logical operation on a 
qubit system, with a given accuracy and in a given time, has been recently postulated. This paper 
is an attempt to formalize this bound and explore the conditions under which it may be expected 
to hold. This is a work in progress and any contributions will be appreciated. 
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It has become of interest lately to explore the con- 
straints that the quantum nature of the control degrees 
of freedom might impose on the practical operation of 
quantum logical gates 0, ||, 0, 1). A very general result 
derived recently by Ozawa ||j is that any quantum gate 
that changes the energy or angular momentum state of a 
qubit will require a minumum number of ancillary bosons 
of the order of 1/e, if it is to have a failure probability 
smaller than e. If the bosons are excitations of a quantum 
harmonic oscillator (such as, e.g., photons) of frequency 
Lo, this becomes a minimum energy requirement 
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in agreement with previous studies ^ which focused on 
the effect of the quantum nature of the electromagnetic 
field on the performance of logical gates. 

Ozawa's result has very wide applicability, but it must 
be kept in mind that it is relatively straightforward (and 
it may be, in fact, advantageous for other practical rea- 
sons) to encode a logical qubit in degenerate states of 
systems of a few qubits, which are mutually intercon- 
vertible without any energy or angular momentum cost: 
for instance, the encoding in a 3-qubit decoherence-free 
subsystem [B uses as the logical zero the state |0)l = 
2-i/2(|01) -niO))|0) of three physical qubits, and as the 
logical one the state \1)l = 6-i/2(|100) + |010) - 2|001)). 
These two states have the same quantum numbers for 
total angular momentum and energy; in fact, they sim- 
ply represent the two different ways to get a state with 
1 = 1/2 and m = —1/2 in a system of three spin-1/2 par- 
ticles. For such an encoding, conservation of total energy 
or angular momentum alone does not appear to restrict 
the possible logical operations. 

I have recently shown Q that in many cases, regardless 
of whether a conservation law is broken or not by the ac- 
tion of the logical gate, there is a minimum requirement 
on the energy of the "control" system, or degree of free- 
dom, of the form (1) if the system is an oscillator, or 
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if the gate is to be carried out in a time T with failure 
probability less than e. My analysis covers gates medi- 
ated by external electromagnetic fields, or by controlled 
collisions between particles, assuming that the fields or 
particles are in minimum uncertainty "coherent states." 
There are, nonetheless, some questions still open, regard- 
ing the full generality of the result, and, for instance, 
whether placing the control degree of freedom in a non- 
classical state (such as a squeezed state) might lower the 
bounds or not. In this paper, which is offered to the 
community as a working document, I shall attempt to 
express the constraint (2) as a formal postulate, and ex- 
hibit a number of worked out examples and ideas for how 
a general proof might proceed. Any help or suggestions 
will be greatly appreciated. 

The reader with only a casual interest may start by 
skipping the detailed calculations of Sections 5 and 6 and 
reading only the more heuristic arguments in the other 
sections. 



II. A FORMAL POSTULATE 

In order to focus only on the constraints arising from 
the quantum nature of the control, and not on those im- 
posed by conservation laws, I have focused on a partic- 
ular kind of two-qubit gate which preserves the qubits' 
energy and angular momentum (assuming the |0) and |1) 
are eigenstates of these variables) , namely, the controlled 
sign-flip gate, which leaves the states |00), |01) and |10) 
unchanged but turns |11) into — 111). The role of the 
control system is, essentially, to switch "on" and "off" a 
Hamiltonian which accomplishes this in a time T, to an 
accuracy given by e. 

To that end, let the control degree of freedom be ini- 
tially in the state and let its self-Hamiltonian be 
Hq. Let the interaction Hamiltonian have the simple 
form Hi = This is the minimal form needed 
for the purpose at hand (question: would it be worth it 
to consider more complicated couplings?) and V need 
depend only on "control" operators. Further, suppose V 
is time-independent in the Schrodinger picture, although 
this may not be necessary. What is necessary is that 
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the interaction be turned on and off only by the con- 
trol system, acting under the influence of its own self- 
Hamiltonian (it follows that |^o) cannot be a stationary 
state). Formally, we require 

(3) 



at the initial and final times, t = Q,T. 



To capture the desired change in sign, we define the 
"failure probability" of the gate by considering what it 
does to a state such as |00) + |11). What we want is 
something like 



i=(|00) + |ll))|^o> - ^(100) - |ll))e-^.C^o*'|^^) (4) 



What we will get is, instead. 



-i=(|00) + |ll))|^o) - i=|00)e-Ho^«orf*'|^^) + _L|ii)e-i Jo"(«o+v)<it'|^^) (5) 



and the "failure probability" can therefore be defined as 
1 minus the square of the overlap between (4) and (5), 
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1 - {,l;o\ei^o"odt'^-^S^(Ho+v)dt'^^^.^ 
1 - (^o|Te-^^o''Vi(t')<it'|^p 
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where the last equation is written in the interaction pic- 
ture, and time-ordering is denoted by T. Now we have all 
the ingredients needed to make a general, formal claim: 
in order to be able to turn on and off an interaction strong 
enough to flip the sign of the wavefunction in (6) over the 
time T, and to do this accurately enough, so that p < e 
(where e is some acceptable error) the state must 
have a minimum energy of the order of 
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This claim involves only the (arbitrary) control system, 
its self-Hamiltonian, and the interaction V. The question 
is, how generally can this be established? 



III. A COUNTEREXAMPLE TO SHOW THAT 
THE CONDITION (3) IS NECESSARY 

At the risk of belaboring the obvious, here is a simple 
example showing that the constraint disappears if one 
allows the interaction to be always "on," that is, if (3) 
does not hold. Let the control degree of freedom be a 
harmonic oscillator, let \ipo) — \n), an energy eigenstate, 
and let V — figa)a. Then one only has to choose T = 
'K / gn and equation (6) will be satisfied exactly, with p = 
0, which means one could make e arbitrarily small, and 
(7) would be violated. 



IV. SOME INSIGHTS ON THE REASON FOR 
THE CONSTRAINT 

At this point it is natural to ask, why, then, is it impos- 
sible to make p in (6) exactly zero if eqs. (3) are enforced? 
A first answer is that, clearly, if V/(i)|V'o) was nearly 
equal to zero at all times the interaction would have es- 
sentially no effect; so then j-i/jo) cannot be an eigenstate of 
Vi{t), which means that Vi{t) will not be sharply defined 
in the state IV'o)- There will be fluctuations, which one 
could formally separate out as 



Vi{t) = {Vi{t)) + /\Vi{t) 



(8) 



with {Vi{t)) = (■(/'ol Vf(t)|V'o)7 and \'4)o) not an eigenstate 
of Ay. In that case, if one chooses V, IV'o), and T so 
that 



{Vi{t'))dt' = irh 



one may estimate the failure probability p given by 
by expanding the exponential, as 
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dt'{M^Vi{t)^Vi{t')\'il;o) (10) 



and this is the approach that was adopted in and in 
most of the examples to follow. (The time ordering has 
been dropped in reaching (|l^), which is presumably not 
important if all one wants is an estimate.) 

On the other hand, perhaps a more fundamental way to 
look at this is to realize that the presence of the operators 
AV^ in the exponent of (|^) changes the state |'0o), so that 
it no longer exactly overlaps with itself, and (|lO|) simply 
estimates the extent of this mismatch. From this point of 
view, what we have here is just the old idea that, in order 
to be able to observe interference in a quantum mechan- 
ical system that is interacting with a classical "appara- 
tus," the apparatus (in this case, the "control" system 



3 



described by \tpo)) must be large enough for the "back 
reaction" of the quantum system on it to be neghgible. 

One way to think along these lines may be as follows. 
Putting together (^) and , one can say that the action 
of the self-Hamiltonian Hq on jV'o) changes it, in a time 
of the order of T, from a state for which V\4>o) ~ to a 
state for which V\4>q) is of the order of {7rfi/T)\tpo) (this 
is a measure of the degree of noncommutativity between 
Hq and V). One may expect the V in the exponent of 
to have a similar effect, then, and "displace" the en- 
ergies of the states making up \ipo) by an approximate 
amount AE ~ ttU/T. The condition (0), then, would 
express the minimum energy that |V'o) must have in or- 
der to still overlap with itself to the degree given by e, 
after its component energies have been "messed up" by 
an amount of the order of AE; in this language, it simply 
reads AE/E < e. 

The next two sections contain a number of worked out 
examples, essentially applying Eq. (|lO|) to various situa- 
tions and showing how something like (0) arises in every 
case. After that, the last Section reexamines the case 
of material particles (where the "control" degree of free- 
dom is, say, a particle's center-of-mass coordinate) in a 
heuristic manner, and presents some final thoughts. 



V. SWITCHING BY MEANS OF A 
(QUANTIZED) E.M. FIELD 

A. Linear coupling 

Assume linear coupling, and a multimode coherent 
state: 
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(12) 



(13) 



Error arises from quantum fluctuations in the ampli- 
tude of the field modes. To keep it smaller than e one 
will require 
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< e 



However, from (|T^) we get that 
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Then (Rj) is only possible if 
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The pulse's "average frequency" is 



and the total pulse energy is 

E field = ^ ^Wfe 



ak\ 



so 
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E field > — 
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Note that for a "static" field, switched on and off over a 
time T, (w) ~ l/T 



B. Nonlinear coupling 

Assume the Hamiltonian is of the form: 

Hi ^hgEP{t)\ll){ll\ (20) 
with g a time- independent coupling constant. Let 
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It is understood that the sum over frequencies is limited 
by the natural frequency response of the system. Let 
£{t) = {i]AE\il;) andE = £ + AE. Then 

AH = hgp£P-^AE 
= ngp£P'' 



k 

The error in the operation of the gate can be estimated 
as 



whereas, on the other hand, we want 
g J £Pdt = TT 



(23) 
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and the left-hand side of this expression can be written 
as 
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£P-^ake-^'^''*dt + c.c. (25) 
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The two conditions 
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and 
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< 3 (27) 



are formally equivalent to ( |13D and (|lj) with a time- 
dependent gk and a modified e, and so the same logic 
applies to yield for the total field energy 



infield > 



(28) 



As long as > 1 (i.e., the Hamiltonian is an analytic 
function of the field) this condition is at least as restric- 
tive as (ITsI). 



C. Squeezing? 

In a coherent state, both "quadratures" of the field- 
amplitude operator ak have the same noise. One could 
imagine a Hamiltonian that couples only to one quadra- 
ture, which could then be squeezed. 

What might happen then could be roughly as follows. 
The fluctuations (squared) in the squeezed quadrature 
would be reduced by a factor e^^"", where r is the squeez- 
ing parameter. This could amount to formally increasing 
e in Eq. (|l4) by a factor e^^. Note that the number of 
photons in the field is now given by |afcp-|-e^'', so, in fact, 
the equation for the field energy might end up reading 



Efield > H^) 
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When this expression is minimized over r, one obtains 
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Note, however, that to couple to a squeezed field one 
typically needs a local oscillator at the carrier frequency 
Lo. Presumably, if w is not sufficiently sharply defined, 
errors in the gate operation will result. This means that 
broadening of lo due to the finite pulse duration must be 
prevented. If the condition 
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is applied to equation (pO[), one obtains again 



E 



field 



> 



eT 



(31) 



(32) 



Thus, it seems that even using squeezing one is still con- 
strained by the inequality (^2|). A more careful study of 
this possibility, however, may be necessary, ideally in the 
context of a specific model for the coupling interaction. 



VI. SWITCHING USING COLLISIONS 
BETWEEN WAVEPACKETS 

A. "Free" particles 

Suppose one arranges to have a collision between the 
two particles involved in the gate operation, with the idea 
that their mutual interaction, V{\ri — r2|), will provide 
the desired phase shift. Work in the center of mass frame 
assuming identical particles; neglect deviations of the 
particles' motion from straight lines at constant speed; 
let b be the distance of closest approach and take that to 
be the x direction. Then what we want is 



i / ^ V{V'^vH'^ + b'^)dt = 7: 



(33) 



where ±v is the y-component of the particles' velocity in 
the CM frame. 

In practice the free wavepackets' x coordinate is un- 
certain by an amount equal to 



Ax{t) = Axo + 



T 
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(34) 



assuming that x and p are initially uncorrelated (at the 
time t = — T/2). This alone causes an uncertainty in (|3^ ) 
whose average magnitude square goes as 



dV dt 
dp p 
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(35) 



with p = (4u^t^+6^)^/-^, and making use of the symmetry 
of the integrands. 

One can now minimize 5"^ with the constraint 
AxoApo > ^/2 (i.e., pick an optimal wavepacket), with 
the resuh Axg = Tn/Am, Ap^ = 2mh/T. Then the 
condition < e becomes 
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Now use (p3|) to eliminate the first factor of lA;, and 



consider the derivative of the left-hand side of (33) with 
respect to the impact parameter b. One easily obtains 
the constraint 
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Vi^/AvH"^ + b'^)dt 



<6 (37) 



This can be simplified a bit by, first, introducing the obvi- 
ous change of variable 2vt — y, with limits of integration 
±?/o — ivT, and then assuming that i/q is large enough 
that no substantial error is introduced by extending the 
integration to ±00. This is reasonable, since the particles 
have to start and end far enough away from each other 
for their mutual interaction to be negligible. We then 
have 
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Viy/y^ + b^)dy 



< e 



(38) 
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The left-hand side of ([S^) is easily evaluated when V{p) 
is any power law, since one can write (j/^ + 5^)~"/^ = 
^ "((y/^)^ + l) "^^ and then the change of variable y/ 6 — 
u results in a factor of 6""-+-'^ times an integral which is 
independent of b. Hence for any n > 1, we get 



2m 52 



< e 



(39) 



and, since, as argued above, we must have b < yo = vT, 
this yields immediately 



mv T 



< € 



or 



mv > 



eT 



(40) 



(41) 



where mv^ is the initial kinetic energy of the two parti- 
cles. 



B. Particles in a harmonic potential 



the following scenario: at time t = we create the two 
wavepackets, a distance AA + b apart, let them oscillate 
towards each other with amplitudes A, so that, at the 
time t — tt/uj, they are closest, a distance b apart; then 
they swing back to the starting position by the time t — 
2'k/lj. At a minimum, one needs to put in and remove 
enough energy to start and stop this pendulum motion. 
Just how this is done is left vague for the moment, but 
it might be important later on. 

In any case, assume that we have an interaction energy 
V{p), as before, but now only in one dimension, with 



(42) 



xi = — I ^ + 2 1 ^ Acosut 



X2 = \ A+ -] + A cos Lot 



and p — X2 — X1 — 2A + b + 2AcosLL}t. The desired action 
is 



1 



2it/u 



V{p)dt = TT 



(43) 



To prevent spreading of the wavepackets during the 
interaction, one could imagine confining the particles in 
a static potential (a time-dependent potential implies a 
time-dependent field, and we are back to the previous 
Section). Assume the potential is harmonic, and consider 



The error operator Ax for a harmonic oscillator is 

A.T = Axq cos cut + Ap sin ujt (44) 
so we find, proceeding as before. 



cos ujtdt I A xo 



dV 
dp 



— — cos Ljtdt ] A xq 
dp 



dV 
dp 



sin ujtdt A po 



(45) 



The contribution of Apo vanishes due to the symmetry of 
the integration, which raises again the possibility of us- 
ing squeezing to improve on the constraint to be derived 
presently. For the moment, however, assume simply that 
a coherent state wavepacket is excited; then Axq is just 
the ordinary zero-point fluctuation of the ground state 
of a harmonic oscillator, A'^xq = h/2muj, and now the 
constraint we have is 




dp 



■ cos ujt dt 



< e 



(46) 



V{p) dt 



The integrals in this case do not seem so easy to evaluate 
in general, but specific cases can readily be done. For 



instance, for a dipole-dipole p ^ interaction one flnds to 
leading order in b 



mcu V 2b , 



(47) 



Note that the energy of each oscillator is ^muj'^A'^, that 
the interaction time T = 2tt/u;, and that, as before, we'll 
want A > 6, so again we find that 



muj'^A'^ > 



eT 



(48) 



As mentioned above, it looks as if one could use a 
squeezed state (squeezed in the position variable) to im- 
prove on the constraint (^). This is because A^pq does 
not appear in Eq. (^5|), which in turn follows from the 
symmetry of the integral over the unperturbed trajectory. 
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Numerical calculations done for a classical particle, how- 
ever, show that (as is only to be expected), as a result of 
the interaction, the particle does not return exactly to the 
starting point. The importance of this mismatch between 
the perturbed and unperturbed wavepackets would only 
be magnified if the quantum wavepacket was squeezed 
in position. Hence, there has to be a limitation to how 
much one can squeeze the position, but it does not seem 
a simple matter to derive it. Specifically, it seems that, 
in the formalism used here, these effects would appear to 
a higher order in the expansion (p5[). 

On the other hand, there may be good (self- 
consistency) reasons to go to higher orders in ( ^5| ) in the 
case of large squeezing. If one has a state that is squeezed 
in position, say, enough to change the dependence of the 
minimum energy on e, from to e~^/^ (the best achiev- 
able in any case, by the arguments of Section V.C), the 
squeezing factor e"^'' in A^xq would have to be of the 
order of -^/e. In that case, the corresponding factor e^'' 
in A^po would be of the order of 1/y/e, and in A^po it 
would be of the order of 1/e. This suggests that in case 
of such extreme squeezing, one would not be justified to 
neglect the higher order terms (in particular, terms of 
order A'*po) in (0)- I currently looking into this. 



VII. DISCUSSION 

For material particles, there is actually a very easy 
way to "derive" the constraint heuristically, based on 
some of the ideas presented in Section 4. The potential 
V produces a force dV/dx on the particle, which, acting 
over a time T, results in a position change (relative to 
the unperturbed wavepacket) of 



dx 



1 dv 
2m dx 



rj-i2 



and a momentum change 



ax 



(49) 



(50) 



From the condition ^ one can estimate V as irh/T and 
dV/dx as Trh/LT, where L is a characteristic length, that 



the particle traverses in the time T (so the velocity v ~ 
L/T). Presumably, the position and momentum change 
will lead to a "mis-overlap" with the original wavepacket 
of the order of {Sx/Ax)'^ and {Sp/Ap)'^, where Ax and 
Ap are the original, intrinsic position and momentum 
uncertainty. The constraint 
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then becomes 
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< e 
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< e 



(51) 



(52) 



\Am? A'^x A^p 
Using the fact that AxAp > h/2 to optimize (minimize 
the left-hand side of (|2|), we find that it reduces to 



2TT^flT 



< e 



which is to say 



1 p n'^h 
—mv > 

2 T 



(53) 



(54) 



if V - L/T. 

The simplicity of this argument is in sharp contrast 
with the complexity of the specific examples worked out 
in the previous section. One must, therefore, ask, is there 
a formally simple way to make the above heuristic argu- 
ment rigorous and general? 

It is probably not hard to see in the case of the electro- 
magnetic field (Section 4) a generalization of the above 
argument to a situation dealing with many "generalized 
coordinates and momenta." One may also ask, at this 
point, whether it is essential, for an inequality of the 
form (2) to hold, that the control system's self-energy be 
a quadratic function of the generalized coordinates and 
momenta. 

Finally, in the case of an electromagnetic pulse the 
pulse duration comes out naturally from the formalism 
of Section 5, whereas for the material particles of Section 
6 there is some ambiguity as to what to use for T. It 
might be nice (necessary?) to define T formally in some 
way, perhaps using the Fourier transform of {Vi{t)). 
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